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Abstract 

The  objective  of  the  proposed  work  is  to  study  the  fundamental  mechanisms  of  active  aeroelastic  tailoring  of  high- 
aspect-ratio  composite  wings.  The  emphasis  will  be  on  actively  exploring  structural  flexibility  to  enhance  flight 
performance  and  reduce  structural  weight,  while  controlling  aeroelastic  instabilities.  To  investigate  the  effects  of 
distributed  anisotropic  strain  actuation  and  their  synergism  with  the  passive  aeroelastic  tailored  composite  structure, 
a  suitable  nonlinear  active  aeroelastic  framework  is  required.  Therefore,  an  analysis  and  design  tool  is  being 
developed  for  high  aspect  ratio  active  flexible  wings.  The  formulation  is  capable  of  modeling  the  nonlinear  large 
deflection  behavior  of  the  wings,  deformation  due  to  applied  voltages  on  the  active  plies  in  the  skin  of  the  wing,  and 
the  unsteady  subsonic  aerodynamic  forces  acting  on  the  wing.  Because  the  wing  is  long  and  slender,  it  can  be 
reasonably  well  modeled  as  a  beam  undergoing  three  dimensional  displacements  and  rotations.  The  cross  sectional 
stiffness,  inertia,  and  actuation  properties  of  the  wing  are  calculated  along  the  span,  and  then  a  1-D  nonlinear  beam 
model  is  constructed.  This  is  a  novel  beam  representation  that  is  based  entirely  on  the  beam  curvatures,  resulting  in 
a  computationally  efficient  low-order  model  suitable  for  preliminary  strucmral  design  and  control  synthesis. 
Preliminary  numerical  results  are  presented  in  this  report. 


ACTIVE  AEROELASTIC  TAILORING  OF  HIGH  ASPECT  RATIO  COMPOSITE 

WINGS 


Abstract 

The  objective  of  the  proposed  work  is  to  study  the  fundamental  mechanisms  of  active  aeroelastic  tailoring  of 
high-aspect-ratio  composite  wings.  The  emphasis  will  be  on  actively  exploring  structural  flexibility  to  enhance  flight 
performance  and  reduce  structural  weight,  while  controlling  aeroelastic  instabilities.  To  investigate  the  effects  of 
distributed  anisotropic  strain  actuation  and  their  synergism  with  the  passive  aeroelastic  tailored  composite  structoe, 
a  suitable  nonlinear  active  aeroelastic  framework  is  required.  Therefore,  an  analysis  and  design  tool  is  being 
developed  for  high  aspect  ratio  active  flexible  wings.  The  formulation  is  capable  of  modeling  the  nonlinear  large 
deflection  behavior  of  the  wings,  deformation  due  to  applied  voltages  on  the  active  plies  in  the  skin  of  the  wing,  and 
the  unstea4y  subsonic  aerodynamic  forces  acting  on  the  wing.  Because  the  wing  is  long  and  slender,  it  can  be 
reasonably  well  modeled  as  a  beam  undergoing  three  dimensional  displacements  and  rotations.  The  cross  sectional 
stiffness,  inertia,  and  actuation  properties  of  the  wing  are  calculated  along  the  span,  and  then  a  1-D  nonlinear  beam 
model  is  constructed.  This  is  a  novel  beam  representation  that  is  based  entirely  on  the  beam  curvatures,  resulting  in 
a  computationally  efficient  low-order  model  suitable  for  preliminary  structural  design  and  control  synthesis. 
Preliminary  numerical  results  are  presented  in  this  report. 

Status  of  Effort 

During  these  first  eight  months  of  this  study,  considerable  progress  was  made  towards  tae  development  of  the 
active-flexible-wing  analysis  formulation.  This  is  a  fundamental  component  for  the  investigation  of  the  effects  of 
embedded  anisotropic  actuation  on  the  aeroelastic  response  of  high-aspect-ratio  composite  wings.  Moreover,  it  was 
found  that  the  nonlinear  beam  problem  can  be  represented  solely  by  tracking  beam  curvatures  in  a  very  natural  way 
without  ever  having  to  introduce  displacement  or  rotation  variables.  These  can  be  later  recovered  by  appropriate 
relations  between  them  and  the  determined  curvatures.  Preliminary  validation  of  the  formulation  has  been 
conducted,  and  further  numerical  studies  are  expected  for  this  next  phase  of  the  program. 

Introduction 

In  the  design  of  aircraft  wings,  the  goal  is  usually  to  minimize  the  weight  of  the  wings  for  a  given  payload 
weight.  In  the  optimal  design  process,  the  structure  of  the  wings  is  reduced  until  design  constraints  become  active. 
Reducing  t^-e  weight  of  the  structure  tends  to  produce  wings  that  are  more  flexible,  resulting  in  a  lower  flutter-speed. 
At  the  (.ptimal  design  point,  it  is  possible  that  one  or  even  all  the  constraints  associated  with  flutter,  divergence,  and 
strength  be  active.  Furthermore,  flexible  wings  will  undergo  large  deflections,  making  their  aeroelastic  response 
more  complex  to  analyze  and  control  using  linear  analysis  tools  and  control  methods.  However,  replacing  the 
conventional  flaps  and  ailerons  with  conformable  surfaces  with  embedded  anisotropic  piezoelectric  materials  in  the 
skin  may  lead  to  lighter  designs.  Two  are  the  main  reasons  for  that: 

1 .  The  actuatuors  are  integral  to  the  strucmre,  and  contribute  to  its  overall  stiffness;  it  eliminates  the  extra  weight 
associated  with,  hydraulic  lines  and  mechanical  linkages. 

2.  The  active  wings  will  be  able  to  extend  the  flutter  and  divergence  boundaries  through  active  control,  allowing  for 
a  new  optimal  design  point  to  be  reached. 

In  the  structural  design  of  active  flexible  wings,  a  daunting  amount  of  trade  studies  must  take  place  in  order  to 
realize  the  best  active/passive  layup  on  the  skin.  The  following  observations  apply: 

1 .  In  order  to  achieve  the  desired  deformation,  the  stiffness  of  the  active  plies  must  be  a  significant  fraction  of  the 
overall  stiffness  of  the  wings. 

2.  Active  ply  orientations  that  produce  high  twist  authority  may  exhibit  other  negative  qualities,  such  as  greater 
static  deflections  and  higher  stresses  in  the  passive  plies. 

3.  The  behavior  of  the  wings  may  depend  highly  on  the  aeroelastic  steady  state  at  the  given  flight  condition. 


To  navigate  such  a  rich  structural  design  space,  where  the  aeroelastic  behavior  is  not  well  known  requires  a 
priori  understanding  of  the  fundamental  effects  involved  in  the  problem  and  potentially  large  degree  of  iteration. 
The  current  phase  of  this  study  concentrates  on  the  development  of  a  numerical  procedure  for  determining  the  best 
active-passive  layup,  sensor  locations,  and  controller  architecture,  which  satisfy  the  lift  and  roll  maneuver 
requirements  over  the  operating  flight  envelope.  The  derivation  of  the  equations  of  motion  for  the  two-wing  active 
roll  model  will  be  briefly  discussed  in  this  report.  Work  on  the  controller/structure  automated  design  scheme  will 
follow. 

Model  Formulation 

The  formulation  of  the  equations  of  motion  governing  the  roll  and  flexible  wing  motion  are  derived  using  the 
principle  of  virtual  work.  The  high  aspect  ratio  wings  are  modeled  as  three-dimensional  beams,  which  can  undergo 
large  deflections.  It  is  required  to  know  the  stiffness,  mass,  and  inertia  distributions  as  function  of  wing  span.  The 
formulation  of  the  equations  of  motion  for  the  wing  consists  of  an  analysis  over  the  cross  section  based  on  the  work 
of  [3].  There,  the  active  and  passive  layup  and  airfoil  shape  yield  a  set  of  actuation,  stif&ess  and  inertia  properties 
for  the  cross  section.  For  a  tapered  wing,  the  cross  sectional  analysis  must  be  performed  at  several  stations  along  the 
span,  such  that  the  properties  can  be  interpolated.  The  result  is  a  one-dimensional  curvilinear  reference  line  with  a 
distribution  of  bending  and  twisting  stiffness,  mass  and  inertia,  and  center  of  gravity  offsets.  The  location  of  the 
beam  reference  line  can  be  arbitrarily  placed  at  a  percentage  of  chord,  or  can  be  placed  at  the  computed  elastic  axis 
or  area  centroid.  The  curvature-based  beam  structural  representation  is  very  efficient  and  appropriate  for  the 
aeroelastic  formulation  being  pursued.  To  the  best  knowledge  of  the  PI,  this  is  a  new  way  of  representing  the  linear 
and  nonlinear  beam  dynamics  problem. 

Kinematics 


Let  =r  [u^  represent  the  body  reference  frame.  The  positions  of  the  points  at  the  root  of  the  right  and 
left  wings  are  fixed  in  the  body  frame  and  are  given  by  £''(0)  and  £^(0),  corresponding  to  the  wing  curvilinear 
coordinate  5=0.  The  direction  vector  ^^(5)  points  along  the  right  wing  reference  axis.  Uy(s)  points  toward  the  tip 
chord,  and  u[{s)  points  in  the  airfoil  z-direction.  The  left  wing  is  similarly  defined.  Then,  assuming  that  the  wing 
reference  axis  remains  of  constant  length,  and  that  shearing  deformation  has  negligible  effect  on  the  kinematics  of 
deformation,  the  positions  and  direction  vectors  as  a  function  of  the  curvilinear  coordinate,  5,  are 


(1) 


where  x^(5),k'^(5),  andAC'j(5)  represent  the  twisting,  flat-wise  bending,  and  chord-wise  bending  curvatures, 
respectively.  The  position  and  direction  vectors  along  the  span  are  given  by 
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Given  the  position  and  direction  vectors  at  the  root,  and  the  distribution  of  curvature  along  s,  the  shape  of  the 
wing  is  completely  defined  through  equation  2.  The  wing  is  discretized  at  nk  equally  spaced  points  at  the 
coordinates  5^ ,  between  which  the  curvatures  are  linearly  interpolated. 


Inertial  Elastic  and  Piezoelectric  Forces  and  Moments 


The  internal  elastic  and  piezoelectric  forces  and  moments  acting  on  a  point  on  the  reference  line  are  given  by 
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where  the  coefficient  matrices  are  found  by  the  cross  sectional  ani’.Iysis  of  [3].  The  4x4  stiffness  matrix  in  this 
derivation  is  a  reduced  form  of  6  x  6  matrix,  in  which  the  transverse  shear  strains  have  been  condensed  out,  but  not 
neglected.  The  formulation  of  [3]  accounts  for  the  effects  of  J»:iiisotropic  piezoelectric  strain  actuation  in  the  active 
composite  plies.  In  the  present  implementation,  the  derivation  of  [3]  has  been  modified  to  facilitate  the  inclusion 
into  a  global  optimization  scheme,  whereby  the  active  and  passive  layup  material  constants,  thicknesses,  ply  angles, 
and  geometry  are  treated  as  variables.  The  number  of  applied  voltages  in  equation  (3)  is  equal  to  the  number  of 
active  plies  in  the  cross  section,  and  will  be  determined  in  the  (semi)  automated  design  process.  The  minor  changes 
to  the  implementation  include  the  data  format  for  specifying  the  cross  section  geometry,  and  the  method  in  which 
path  integrals  are  carried  out.  The  present  implementation  was  verified  successfully  against  results  of  [3]. 


During  the  cross-section  calculations,  a  set  of  recovery  matrices  are  constructed  which  relate  the  1-D  beam  strains 
and  applied  voltages  to  the  distribution  of  strain  and  stress  in  each  ply  of  the  cross  section.  These  values  will  be  an 
important  driving  force  in  the  design  optimization  to  influence  changes  in  the  cross  section  variables.  The  ply 
strains  are  recovered  at  finely  discretized  points  on  the  cross  section  contour  and  the  ply  stresses  are  determined 
based  on  the  plane  stress  relation, 
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where  the  free  strain  due  to  actuation  is  taken  into  account  in  the  active  plies.  The  present  implementation  allows 
for  the  inclusion  of  an  internal  web  which  may  be  active. 


Inertial  Forces  and  Moments 


The  inertial  forces  and  moments  acting  at  the  center  of  gravity  of  the  cross  section  must  be  referred  to  the  beam 
reference  line.  The  acceleration  of  the  center  of  gravity  of  the  cross  section  is  given  by  ==  ^  where  the 

subscripts  eg  and  ra  correspond  to  the  center  of  gravity  and  reference  axis,  respectively. 

Then  the  forces  and  moments  per  unit  length  acting  on  the  reference  line  are  related  to  the  reference  line 
accelerations  by  the  matrix  relationship 
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where  the  forces  and  accelerations  are  resolved  in  the  local  frame  of  the  deformed  reference  line.  A  vector  quantity 
can  be  converted  from  the  local  to  body  frames  by  the  rotation  matrix  f  =  It  is  more  convenient,  however, 

to  construct  the  matrix  that  converts  all  of  the  nodal  force  vectors  from  the  local  to  the  body  frame  at  once.  Then  a 
vector  of  3np  x  1  nodal  forces  is  converted  from  the  set  of  np  local  frames  to  the  body  frame  by  the  relation 
=  where  is  a  2  npx  3np  rotation  matrix. 


Equations  of  Motion 


To  apply  the  Principle  of  Virtual  Work,  the  virtual  work  done  by  all  internal  and  external  forces  must  be  expressed 
in  terms  of  the  set  of  independent  generalized  variables.  The  internal  forces  include  elastic  and  actuation  forces 
acting  on  the  reference  line.  The  external  forces  are  all  referred  to  the  reference  line  and  include  gravity,  inertial, 
and  aerodynamic  forces.  The  total  work  done  by  these  forces  can  be  divided  into  three  parts, 

1.  The  virtual  work  done  on  the  right  wing: 

=  J{5£'(F“”  -[/:]£) +  5P'(£‘^  -[p]£-[prl)  +  50'(M*  +  [pr]£-[/l}rf5 

right  wing 

2.  The  virtual  work  done  on  the  left  wing: 

5W^  =  -[KW  +  SP'iF^' -[p]E-lpr]e)  +  S9'iM'^+[pr]P-[I]e]ds 

left  wing 

3.  The  work  done  against  the  fuselage  inertia  (roll  only) 

5^3  =  S<l>J^(p  (8) 

where  (p  is  the  roll  angle,  and  is  the  inertia  of  the  fuselage.  The  local  accelerations  are  related  to  the  generalized 

accelerations  through  the  Jacobian  matrix  relationship, 

tu..,  =  +  [^]?>)  (9) 

Discretization  and  Numerical  Itegration 

The  curvature  variables  are  discretized  at  nk  equally  spaced  points  along  the  wing,  at  the  coordinates,  The 
displacements,  rotations,  and  forces  and  moments  are  discretized  at  np  equally  spaced  points  along  the  wing,  at  the 
coordinates,  5^.  Since  a  linear  distribution  of  curvature  cannot  be  represented  by  a  linear  distribution  of 

distplacements,  the  displacement  discretization  is  chosen  to  be  finer  than  the  curvature  discretization. 

A  vector  of  forces  defined  at  the  coordinates  can  be  interpolated  to  the  coordinates  5y-by  the  matrix  relation 
s^-[ap]Sj,^  The  coordinates  5^ are  of  fine  spacing  suitable  for  numerical  integration.  The  integral  equations 

describing  virtual  work  can  be  written  as  matrix  relations  by  using  interpolation  matrices  and  inner  products.  In 
matrix  form,  the  total  virtual  work  relations  becomes  - 

1.  Right  wing 
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2.  Left  wing 
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Application  of  the  PVW  yields  the  second  order  matrix  equilibrium  equations  of  motion,  which  includes  the 
flexible  motion,  roll  motion  input  at  the  wing  root,  and  applied  voltages,  forces,  and  moments. 
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The  aerodynamic  and  gravity  forces  and  moments  expressed  in  the  input  vectors  on  the  right  hand  side  must  be 
expressed  in  terms  of  the  generalized  variables  to  produce  an  effective  state-space  model  of  the  system. 

Unsteady  Aerodynamics  Model 


The  unsteady  aerodynamic  forces  and  moments  are  computed  using  the  approach  of  [5],  which  calculates  the 
unsteady  loads  on  a  deformable  airfoil  in  subsonic  flow.  The  aerodynamic  loads  can  be  represented  in  terms  of  the 
well-known  pitch  and  plunge  variables,  as  well  as  the  induced  flow  due  to  free  vorticity.  Tlie  method  allows  for  the 
aerodynamic  loads  to  be  efficiently  incorporated  into  the  state  space  representation,  and  requires  induced  flow  states 
to  be  added  to  the  structural  states.  By  relating  the  pitch  and  plunge  motion  to  the  generalized  variables,  the 
aerodynamic  forces  can  be  written  in  matrix  form  as 

Mr"'  =  ^nlLr  +  ^13^^  +  ^lA^r  +  ^16^  +  ^110  Ar  + 

M""'  =  ^22^  +  ^23^^  +  +  ^26^  +  ^211^/  +  ^213& 

where  A  and  a  are  die  inflow  and  angle  of  attack  vectors  defined  at  the  coordinates,  5^.  Mr"'»£r"'» 

Mr""  are  vectors  of  aerodynamic  forces  and  moments  on  the  right  and  left  wings,  respectively.  The  method  can  be 

extended  to  include  dynamic  stall  effects  using  the  ONERA  approach. 

Aeroelastic  Model 

By  combining  the  aerodynamic  forces  with  the  structure  model,  the  full  set  of  equations  describing  the 
aeroelastic  motion  is  constructed, 
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where  SKq  Kq  is  the  curvature  vector  of  the  linearized  system,  K;  is  the  initial  curvature  vector  of  the 

wing,  and  5k  is  the  change  in  the  curvature,  is  the  direction  of  gravity  in  the  body  reference  frame.  The  system 
matrices  depend  on  the  current  state  and  the  flight  condition  and  must  be  continuously  updated. 

Preliminary  Numerical  Study 

The  present  study  investigates  the  effects  of  the  active  ply  angles  on  the  static  and  dynamic  behavior  of  a  wind- 
tunnel  size  model  of  an  active  flexible  wing.  The  analysis  model  is  a  NACA0014  airfoil  with  four  (passive)  plies 
(0°/90°)  of  E-glass  running  the  full  length  of  the  chord  on  both  the  top  and  bottom  wing  surfaces.  Also,  tw  o  active 


fiber  composite  (AFC)  plies  at  +9  and  -9  on  each  of  the  top  and  bottom  surfaces,  covering  approximately  95%  of  the 
surface.  The  AFCs  are  used  here  as  a  realization  of  an  anisotropic  piezocomposite  actuator.  The  angles  of  &e 
passive  plies  are  held  constant.  The  wing  has  a  2m  semi-span,  with  a  root  chord  of  15  cm.  The  wing  starts  tapering 
linearly  at  0.5  m  from  the  root,  reaching  a  tip  chord  of  10  cm.  Actuating  the  +9  plies  with  positive  voltage  and  the 

—9  plies  with  negative  voltage  produces  twist.  Bending  actuation  is  achieved  by  applying  positive  voltage  to  the  top 
plies  and  a  negative  voltage  to  the  bottom  plies. 

Figure  2  shows  the  effectiveness  of  the  actuators  for  producing  twist  as  a  function  of  AFC  ply  angle.  A  2000  V 
twist  input  was  applied.  At  0®  the  actuators  have  no  capability  for  producing  twist,  and  the  maximum  effectiveness 
is  shown  to  be  at  45°. 

The  increase  in  twist  capability  comes  at  the  expense  of  bending  authority.  The  0°  plies  are  much  more 
effective  in  producing  bend  deformation  than  the  45°  plies.  As  seen  in  Figure  3,  which  plots  the  maximum  bending 
curvature  for  ply  angles  from  0°  to  50°,  the  0°  plies  produce  more  than  three  times  the  curvature  as  the  45°plies. 

Variations  in  active  ply  angle  that  benefit  twist  authority  may  have  an  adverse  effect  on  the  static  deflection  due 
to  aerodynamic  loading,  which  in  turn  has  an  effect  on  the  natural  frequencies.  Note  that  the  natural  frequencies  of 
the  wing  do  not  only  depend  on  the  mass  and  stiffness  distribution,  but  also  on  the  steady  deformed  shape  due  to 
aerodynamic  loading.  When  the  wing  is  bent,  the  first  torsion  and  first  chordwise  bending  modes  tend  to  coalesce 
into  a  "rocking"  mode.  This  behavior  caimot  be  capmred  using  a  linear  model,  linearized  about  the  undeformed 
state. 

As  seen  in  Figure  4,  which  plots  the  wing  tip  displacement  for  flight  speed  U  =  25  m/s  as  the  angle  of  the  active 
plies  is  varied  from  0°  to  45°  the  tip  displacement  increases  by  22%. 

Figure  5  plots  the  movement  of  the  first  four  natural  frequencies  with  respect  to  the  active  ply  angle.  For  each 
angle,  the  wing  was  brought  to  its  steady  deformed  shape  at  U  =  30  m/s,  and  then  the  natural  frequencies  were 
obtained.  Figure  6  shows  the  reduction  in  flutter  speed  as  the  active  ply  angle  is  varied. 

Those  simple  results  exemplify  the  capabilities  of  the  proposed  analysis  and  showed  that  the  choice  of  actuation 
layup  variables  has  a  profound  effect  on  the  aeroelastic  behavior  of  active  conformable  wings.  By  investigating 
only  changes  in  the  actuation  ply  angles,  it  was  shown  that  adjustments  to  improve  the  twist  authority  of  the  active 
plies  may  lead  to  a  deterioration  in  aeroelastic  performance  for  reasons  that  are  difficult  for  the  designer  to 
understand  without  iterative  data. 

A  redr  otion  in  open  loop  flutter  speed  due  to  increasing  the  active  ply  angles  to  45°  may  be  recovered  when 
active  control  is  employed,  so  the  closed  loop  response  of  the  wing,  as  well  as  static  performance,  must  be  the 
measure  of  the  merit  of  a  particular  wing  structure. 
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Figure  1 :  Two-wing  roll  model  showing  gravity  and  aerodynamic  loading  distributions 


Figure  2:  Changing  in  tip  twist  angle  amplitude  for  different  actuator  ply  angle  in  the  active  wing 

(2000  V  input) 


Figure  3;  Changing  in  maximum  bending  curvature  for  different  actuator  ply  angle  in  the  active 

wing  (2000  V  input) 
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Figure  4:  Changing  in  vertical  tip  displacement  for  different  actuator  ply  angle  in  the  active  wing 

(U=25  m/s) 
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Figure  5;  Variation  in  the  natural  frequencies  of  the  active  wing  due  to  changing  in  the  actuator 

ply  angle 


OPEN  LOOP  FLIJITER  SPEED  V3.  ACTIVE  PLY  ANGLE 


Figure  6:  Reduction  on  the  linearized  flutter  speed  of  the  active  wing  due  to  changing  in  the 

actuator  ply  angle 
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Figure  7:  Wing  deformation  due  to  2000  V  bending  actuation  (figure  compares  deformed  and 

undeformed  states) 
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